General entanglement scaling laws from time evolution 
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We establish a general scaling law for the entanglement of a large class of ground states and dynamically 
evolving states of quantum spin chains: we show that the geometric entropy of a distinguished block saturates, 
and hence follows an entanglement-boundary law. These results apply to any ground state of a gapped model 
resulting from dynamics generated by a local hamiltonian, as well as, dually, to states that are generated via 
a sudden quench of an interaction as recently studied in the case of dynamics of quantum phase transitions. 
We achieve these results by exploiting ideas from quantum information theory and making use of the powerful 
tools provided by Lieb-Robinson bounds. We also show that there exist noncritical fermionic systems and 
equivalent spin chains with rapidly decaying interactions whose geometric entropy scales logarithmically with 
block length. Implications for the classical simulatability are outlined. 

PACS numbers: 



At the heart of the intriguing complexity of describing 
quantum many-body systems is the entanglement contained in 
the system's state: if the state is highly entangled, one needs 
a large number of parameters to describe it classically. The 
scaling of the geometric entropy [1-1 1] - the degree of entan- 
glement of a distinguished subsystem with respect to the rest 
- for quantum many-particle systems, such as those encoun- 
tered in condensed matter physics, is the crucial parameter 
which quantifies whether the state is hard or easy to simulate 
using density-matrix renormalisation group methods |8]. 

Recently, motivated partially by questions of simulatabil- 
ity, there has been a considerable effort to precisely charac- 
terise scaling laws for ground- state entanglement, which we 
call the static geometric entropy [1-11] Indeed, substantial 
progress has been made in answering this difficult question: 
earlier conjectures, for which there was only numerical evi- 
dence, could be resolved. For example, it is now known that 
for gapped bosonic harmonic systems, such as free field mod- 
els III2I1 . the geometric entropy scales like the boundary area 
of a distinguished region, and not the volume The only 
precise results available at the current time pertain to quasi- 
free (or Gaussian) bosonic and fermionic models 0001 and 
equivalent ID spin chains. Apart from integrable systems and 
matrix-product state hamiltonians (which satisfy an area law 
by construction (3), there is a dearth of results concerning 
static geometric entropy for systems as simple as the ID spin- 
1 Heisenberg model. How does the geometric entropy scale 
for general interacting systems? 

There are also very few results available about the strongly 
related case of geometric entropy for dynamically evolving 
states (31. The dynamic geometric entropy occupies centre 
stage when trying to simulate systems which undergo a sud- 
den quench of a local interaction, for example, when a system 
is in a Mott phase when the hopping is suddenly altered. In 
the Mott phase t he g eometric entropy is zero and grows as the 
system evolves . It is far from obvious how the geometric 
entropy should scale as a function of time in these and similar 
systems dynamically undergoing a quantum phase transition. 



In this Letter we establish the first scaling laws for the 
geometric entropy of a general class of quantum states that 
goes significantly beyond Gaussian models. On one hand, we 
will show that if any state of ID spins whose geometric en- 
tropy satisfies a boundary law (i.e., it saturates as a function 
of n, the number of spins) is subjected to dynamics accord- 
ing to an arbitrary ID local model K for any constant time t 
then the dynamic geometric entropy will continue to satisfy a 
boundary law, albeit saturating at a larger constant which de- 
pends linearly on |t | . On the other hand, when considering the 
time evolution generated by the local hamiltonian K, the state 
that results from this time evolution can be thought of as the 
ground state of a gapped hamiltonian, local or with rapidly 
decaying interactions f\3\. All constituents will eventually 
become correlated, but the entanglement built up between re- 
mote parts can be bounded, an intuition that we will cast into 
a rigourous form. Hence, this reasoning is a device that allows 
us to establish the result that the static as well as the dynam- 
ical geometric entropy of a large class of models, including 
strongly correlated systems, satisfies a boundary law. 
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FIG. 1: Hierarchy of unitaries VKi (t) , VK2 (t) , ... with exponentially 
decreasing entangling power over the boundary of A and B . 

To actually carry out the argument outlined above we use 
the powerful machinery of Lieb-Robinson bounds fi^ fnlfTsL 
191 . The intuition we develop is that in a many-body system 
with local interactions there is a finite speed of sound, and 
hence a finite velocity of information transfer, resulting from 
local interactions. The Lieb-Robinson bound is the precise 
quantification of this statement: it says that the norm of the 
commutator of two operators, one of which is evolving ac- 



2 



cording to local dynamics, is exponentially small in the sep- 
aration between the two operators for short times. This in- 
equality allows us to precisely bound the entanglement that 
can develop across the boundary of a distinguished region for 
short times. In turn, we find that for large times of the order 
of the logarithm of the number of spins, the boundary law for 
the dynamic geometric entropy breaks down. We show this 
dually by explicitly constructing a local translation-invariant 
gapped system whose ground state violates an area law. 

Geometric entropy in spin chains. - We will, for the sake 
of clarity, describe our results mainly for a finite chain C of 
n distinguishable spin- 1/2 particles. The family H of local 
hamiltonians we focus on (which implicitly depends on n) is 
defined by = XI j=o ' where Hj acts nontrivially only 
on spins j and j + 1. We set the energy scale by assuming that 
\\Hj II scales as a constant with n for all j = 0, 1, . . . , n — 1, 
where || • || denotes the operator norm. The interaction terms 
Hj can, w.l.o.g., be taken to be positive semidefinite, and may 
depend on time as Hj = Hj (t). 

Consider a bi-partition of the chain into two contiguous 
blocks A and B of spins of sizes m = \A\ and n — m = \B\, 
m < n. We will find boundary laws - a saturation of the 
block entanglement - independent of the system size (we 
avoid the technicalities arising in the case of infinite systems 
which might obscure the main point). For simplicity we as- 
sume m < n/2 and we let |V^(t)) = e*^^|0). The initial state 
is taken to be a product state |0), but the argument is general 
enough to be applicable for any matrix-product or finitely cor- 
related 10, li^] initial state, or a state resulting from a quantum 
cellular automata. Consider the Schmidt decomposition 

where the Sj {t) are the non-increasingly ordered Schmidt co- 
efficients. They are given by the eigenvalues of M{t) = 

C{t)CHt), where m)) = E'Zo' Et7~' C,,km,k), 
and the \j) and \k) form an orthonormal basis for Ha and 
Hb, respectively. The geometric entropy of a block A, or 
the block entanglement, is given by the von-Neumann en- 
tropy 6'(m) = — XI^-^q ^3^^^2{^3) 12111 • We denote by 

Ha = Zlj=o^ = Yl^Zm Hj ^^^^^ /?arr^ of 

the hamiltonian H, which act nontrivially only on subsystem 
A and B, whereas Hi = Hm-i denotes the interaction term. 

Entanglement scaling in dynamically evolving quantum 
states. - In this section we prove an upper bound for the dy- 
namic geometric entropy S{m), S(m) < cq + ci\t\ , where 
Co, ci > depend only on \\h\\ |22] and not on n. Thus, the 
entropy of the block A scales, asymptotically, less than a con- 
stant. Our first step is to obtain the decomposition 

e^'"' = {UA{t)^UB{t))V{t). 

We do this by guessing UA{t) = e'*^^ and Usit) = e'^^^ . 
The idea here is that the dynamics generated by H should be 
similar to those generated by Ha + Hb preceded by a uni- 
tary V{t) that "patches up" the removed interaction. We ob- 
tain a differential equation for l/(t) = e-^^(^^+^^)e'^^ = 



^-itiH-Hj)^itH. dV{t)/dt = V{t)L{t),. The "hamilto- 
nian" L{t) = iHi + J^r^{[H,Hi])du, where r^^ {N) = 
^-itM j^^itM ^Qj. operators M, is antihermitian, so that the 
dynamics this integro-differential equation generates is uni- 
tary. 

Our strategy at this point is to decompose V{t) into a prod- 
uct of strictly local unitary operations ll23ll Va- (t) which act 
nontrivially only on Aj = {x : d{x^ m) < j}, which consists 
of only those sites within a distance j from the boundary. This 
decomposition for V{t), depicted in Fig. 1, is then 

Vit) = Wn-mit)Wn-m-lit) • • • W2{t)Wi{t), 

where Wk{t) = Va^ (t)V^^_^ (t) acts nontrivially only on 
A/e. Also, we set Wi{t) VXi(t) and Wn-m{t) 
V{t)V^_^_i{t). Physically, we expect that the unitary op- 
erators Wi {t) are successively weaker and weaker. To find a 
bound on || (t) — I|| we now invoke the machinery of Lieb- 
Robinson bounds 1 16] (see Ref. ifiol] for a simple direct proof) 
on the speed of sound in systems evolving according to local 
dynamics: the strongest available such bound fioll yields 

||rf^HM)-rf"-^(M)||<J,|t|7/!, 

with 5i = ||M||2^||/i||^ where M = [Ha,,Hi] is the same 
for all k as long as /c > 1. This is indeed the powerful tool 
we need to derive the desired bound concerning the deviation 
from each of the unitaries Wi (t) from the identity: it exponen- 
tially bounds the information spread in a system undergoing 
dynamics under a local hamiltonian. We hence find 

\mt) - III = \\V^,{t) - V^,_M\ < Si\trV{l + 2)!. 

In the last inequality, we have expressed the operators Va^ {t) 
as integrals in time t2M . We note that this bound is decaying 
faster than exponential in This bound tells us that we can 
write V^^(t) =I^ei{t)Xi{t), where \\Xi{t)\\ = land£^(t) < 
min{2, Si\t\^~^'^ / {I + 2)!}. Let us now consider the action of 
^ztH initial product state vector |0), 

n—m 
k=l 

We define |t/-;(t)) = Ylk=iWk{t)\0). Let US now choose / 
large enough so that this bound for Wi (t) is strong enough that 
ei{t) < min{2, ||M||(^^|t|^+V(/ + 2)!} is small [25]. Then 
\il^i{t)) is in a product with respect to the spins outside the re- 
gion Ai There are, in general, at most 2^ nonzero Schmidt co- 
efficients for {ipi {t)) with respect to the bi-partition AB. Now 
we consider the action of Wi-^i{t) on \il^i{t)) which yields 
|?/^/+i(t)). In the computational basis, we have 

2^ + 1-1 
j,k=0 

setting IMt)) = T.tlLo^{Ci)j^k{t)\j,k) and 

Xi^,{t)\Mt)) = E-^^o'(A),,fcW|j,^). The nor- 
malisation condition and ||X^(t)|| < 1 imply that ||Q|| < 1 
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and < 1. We use Weyl's perturbation theorem l28l] 

to bound the Schmidt coefficients s^^^^^ of given 
by the eigenvalues of {Q + ei^i{t)Di){Ci + £z+i(^)A)^ 
We apply Weyl's perturbation theorem to the operators 
P = CiCj andQ = ei^i{t){CiDl ^ DiCj) ^ ef^,{t)DiDj 
with HQ II < c£/+i(t), with c > 0. The eigenvalues of P 
are precisely the 2^ Schmidt coefficients of |V^z(t)). Weyl's 
perturbation theorem tells us that the first 2^ eigenvalues 
of P + Q have to be close to the Schmidt coefficients of 
\ipi{t)) and the remaining 2^ eigenvalues have magnitude 
less than ei-^i{t) O- Exploiting these bounds iteratively, 
we find that the Schmidt coefficients satisfy the bound 
Sj{t) < min{l/2'^l^l,2'^l^l"''-^'}, for some > 0. Hence 
the geometric entropy S{m) satisfies the upper bound 
S{m) < Co + ci|t|, where co,ci > 0. This holds true for 
all n. In other words, we can perform "the limit of infinite 
system size" n ^ oo. When we let |t| = log(n) our bounds 
begin to fall apart: the Lieb-Robinson bound becomes a 
polynomial bound. This situation can be saturated, see below. 

Entropy -boundary laws for approximately local quantum 
spin systems. - We now show that the entropy-area law for dy- 
namically evolving product states implies entropy-area laws 
for the ground states of noncritical approximately local quan- 
tum spin systems. The product |0) is the unique ground state 
of the hamiltonian Z = — Y^^Zl ^j- Let H be our hamil- 
tonian. Then \ip{t)) is the unique ground state of the new 
hamiltonian K = e^^^Ze"*^^, having exactly the same spec- 
trum as Z. Moreover, while K is no longer strictly local 
in general, it is approximately local with exponentially de- 
caying interactions. The way to see this is to apply a Lieb- 
Robinson bound to the interaction term e*^^cr|e~*^^: We 
consider the difference between e*^^cr^e~*^^, having sup- 
port equal to C, and the strictly local e'^^^^kU) a^e'^'^^^kU) ^ 
with Ak{j) = {x : d{x^j) < k}, which has support on 
2k -\- 1 sites. This difference can be bounded using the Lieb- 
Robinson bound, Wn^icr^) - rf^'^'^(o-|)|| < ce'^l^l-^^, with 
c > 0. Thus, the interaction term (a^) couples spins 
from site j exponentially weakly. What sort of hamiltonians 

_ ^itH ^^-itH _ clearly a large class of gapped models - 
arise in this way? Insight can be provided by the following 
example: Let H = ajaj^^. For small t the hamiltonian 
K will look like K = Z^\(t) ^j^j^i +0(^2). 
In this case K is similar to the XY model in an external mag- 
netic field with small higher order terms. With the inclusion 
of larger neighborhoods, in turn, local hamiltonians can be 
approximated to any accuracy. Another useful hamiltonian 
which can arise in this way is the strictly local cluster hamilto- 
nian 1 29] (carrying over also to the higher-dimensional case): 
set H = cF^cF^j^-^. In this case, when t = 7r/2, is the 
hamiltonian having the cluster state as a unique ground state. 

Logarithmic divergence of geometric entropy of gapped 
systems. - We now construct an explicit situation where a 
gapped ID spin system C ndeed violates the entanglement- 
boundary law. We again consider a family of spin systems, 
consisting of a block A consisting of m spins, and a block B 
containing the remaining n — m spins. As before, S{m) is 



defined to be the geometric entropy of a block A of m spins 
in the limit of an infinite chain, for simplicity with periodic 
boundary conditions. By virtue of the familiar Jordan- Wigner 
transformation iEtIi . we may consider the fermionic model 

n-l 

where Mi G M, / = 0, . . . , n — 1. The hermiticity of H and 
the periodic boundary conditions are reflected by the condi- 
tions Ml = M-i for afl / = 0, . . . , n - 1 and Mi = Mi^n- 
We can easily map the above hamiltonian onto the one for 
non-interacting fermions, preserving the anti-commutation re- 
lations: H = YlkZo ^kbl.bk, where A: = 0, . . . , n — 1, are 
the eigenvalues of M, given by ek = EjJo e^^^^^+^^^/^M,. 
The ground state can then be easily found: it is the state with 
unit occupancy for each k with Sk < 0. If the value is not 
contained in the spectrum, this ground state is non-degenerate. 
We now consider the subsystem A. The reduced state pm of 
this block is characterised by the spectrum of the real symmet- 
ric m x m Toeplitz matrix 112^ which defines the second 
moments of fermionic operators (41 The l-th row of 

this matrix is given by t_/+2, to, . . . , tm-i), where 

ti = YlkZo e~*^^£/e/(^|£/c|). At this point, we may take the 
limit n ^ oo, for fixed m, and consider long-ranged interac- 
tions, and hence sequences of couplings {Mj^^p^, Mi G M. 
This means that in the continuum limit, we can consider func- 
tions ^ : (0, 27r] R, representing the spectrum of the in- 
teraction matrix, and t^ = 1/{27t) J^^ e~'^^^(f){x) /\(l){x)\dx. 
We can now make use of a very useful bound of Ref. 101, 
stating that S{m) > -(log2 | det[T^]|)/2. Hence, to show 
that S{m) = 17 (log m), we have to bound the Toeplitz 
determinant det[T^]. This we can do usinga proven in- 
stance of the Fisher-Hartwig conjecture (Ifl UM^ determin- 
ing the scaling of the determinants of Toeplitz matrices. Us- 
ing these ideas, we are now construct a model with the men- 
tioned surprising properties: We take the interactions Mk = 
J^"" (l){x)e-'''^/{27T) to be given by Mj, = -z(e^^^/2 _ 
1)^(1 + e*^^/2)/(2e2^^^A;7r), so a 1/k decay of the interac- 
tions, as in case of an unshielded Coulomb interaction. This 
gives rise to the Fourier transform (j) that takes the value 
1 in X e (0,7r/2], and (37r/2,27r] and the value -1 in 
X G (7r/2, 37r/2]. In this setting, the proven instance of the 
Fisher-Hartwig conjecture then indeed allows us to argue that 
|det[T^]| = Q{\ogm) 1I3TI] . This hamiltonian is obviously 
gapped: the quasi-particle excitation spectrum is even con- 
stant, and never crosses zero, so it defines a gapped system. 
Still, we find a logarithmically divergent geometric entropy. 
This is an example of a ground state that is not covered by the 
above statement for small times. 

Outlook. - In this work we have introduced an approach to 
assess geometric entropies in many-body systems. We have 
found that many ground states of quasi-local gapped hamilto- 
nians, while being far from quasi-free, still exhibit a saturat- 
ing geometric entanglement, and hence an entanglement- area 
law. The studied gapped systems are rigorously classically 
efficiently simulatable: one can obtain all expectation values 
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of local observables with polynomial computational resources 
tnll . Simulatability is closely linked with ID entropy bound- 
ary laws 0]. This connection is even more direct in our case 
because matrix-product states which faithfully approximate 
our ground states can be explicitly constructed LI 81 . Such ef- 
ficient descriptions in terms of matrix-product states would 
also be generated by an eventually successful application of 
the DMRG algorithm to our systems. 

Two-dimensional systems are in principle accessible with 
the methods introduced here. This method opens up the way 
toward studying the complexity of gapped many-body sys- 
tems and the accompanying ground- state entanglement scal- 



ing (as well as capacities of quantum channels based on inter- 
acting systems 1331). Intriguingly, we finally found an exam- 
ple of a gapped system with a divergent block entanglement, 
rendering the connection between criticality and validity of an 
area theorem more complex than anticipated. 
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